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In this work we have used 3D hydrodynamical (C0 5 BOLD) and ID hydrostatic (LHD) stellar 
atmosphere models to study the importance of convection and horizontal temperature inhomo- 
geneities in stellar abundance work related to late-type giants. We have found that for a number 
of key elements, such as Na, Mg, Si, Ca, Ti, Fe, Ni, Zn, Ba, Eu, differences in abundances 
predicted by 3D and ID models are typically minor (< 0.1 dex) at solar metallicity. However, 
at [M/H] = -3 they become larger and reach to -0.5 . . . -0.8 dex. In case of neutral atoms and 
fixed metallicity, the largest abundance differences were obtained for the spectral lines with low- 
est excitation potential, while for ionized species the largest 3D-1D abundance differences were 
found for lines of highest excitation potential. The large abundance differences at low metallicity 
are caused by large horizontal temperature fluctuations and lower mean temperature in the outer 
layers of the 3D hydrodynamical model compared with its ID counterpart. 
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1. Introduction 

Red giants are important tracers of intermediate age and old stellar populations and due to their 
high luminosity they are available for spectroscopic study beyond the Milky Way. Their abundances 
are generally determined using stationary ID stellar atmosphere models which treat convection in a 
parametric way (e.g. microturbulence velocity, mixing length). These shortcomings of the classical 
ID models are overcome in 3D hydrodynamic stellar atmosphere models which treat convection 
by solving time dependent hydrodynamic equations. Due to different approaches in the treatment 
of convection one may expect differences in the predicted strengths of spectral lines, and thus - 
differences in elemental abundances derived with the two types of models. 

Comparison of abundance differences arising between the predictions of 3D and ID stellar 
atmosphere models for red giants was made by [4]. Here we extended their work and used the 3D 
hydrodynamic stellar atmosphere models to investigate the influence of convection and horizon- 
tal temperature inhomogeneities on the predicted spectral line strength and thus - on the derived 
elemental abundances in red giant stars. For this purpose we utilized a broader list of chemical 
elements and spectral line parameters, such as wavelength and excitation potential of the lower 
level, %. This analysis was performed by means of the 3D-1D abundance corrections defined as 
an abundance difference between the 3D and ID curves of growth obtained at a given equivalent 
width. The concept of the 3D-1D abundance correction was described in more detail in [3]. 



2. Models 

In this work we used 3D hydrodynamic and ID hydrostatic models calculated using the 
C0 5 BOLD 1 [5,6] and LHD [2] codes, respectively. Both 3D and ID models shared identical atmo- 
spheric parameters (Table 1), chemical composition, equation of state and opacities. Solar scaled 
chemical composition with alpha-element enhancement of +0.4 dex for metallicities [M/H] 2 < -1 
is assumed in both 3D and ID models [7]. Temperature stratifications of 3D and ID models for 
[M/H] = and -3 used in the current study are shown in Fig. 1. 
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Table 1: Atmospheric parameters of the CO s BOLD models used in this work. 



1 http://www.astro.uu.se/~bf/co5bold_main.html 

2 [M/H] = log[N(M)/N(H)]+ - log[N(M)/N(H)] . 
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Figure 1: Temperature profiles in two late-type giant models plotted versus Rosseland optical depth (left: 
^eff/log g/[MM] = 4970/2.5/0.0; right: 5020/2.5/-3.0). The temperature distribution in the 3D model is 
shown with color-coded temperature increasing from blue to red. The dashed line shows the mean temper- 
ature of the 3D model averaged on surfaces of equal optical depth. The solid line is the temperature profile 
of the ID model. Horizontal bars mark the regions where 90% of the equivalent width of the Fe I and Fe II 
spectral line at X = 4000 A and % = eV and 6 eV is forming in the 3D (black) and ID (blue) atmosphere 
models. 



3. Spectral line synthesis 

Spectral line synthesis was done using the LINFOR3D 3 code which solves the 3D radiative 
transfer problem under the assumption of LTE. 

The 3D-1D abundance corrections were calculated for a number of neutral and singly ionized 
elements: Na I, Mg I, Mg II, Si I, Si II, Ca I, Ca II, Ti I, Ti II, Fe I, Fe II, Ni I, Ni II, Zn I, Zn II, 
Ba II and Eu II. 

Only weak lines (W < 5 mA) were used in calculations in order to avoid the influence of 
microturbulence parameter on the strength of ID model lines. To investigate the dependence of 
3D-1D abundance correction on wavelength and excitation potential we have synthesized fictitious 
lines at A = 4000 A and 8500 A and excitation potentials in steps of 2 eV, from to 6 eV. 



4. Results 

4.1 Abundance corrections for neutral atoms 

The dependence of the 3D-1D abundance corrections on metallicity for neutral atoms is shown 
in Fig. 2 and 3, at X = 4000 A and 8500 A, respectively. At solar metallicity the abundance 
corrections are small for all investigated chemical elements, never exceeding ±0.05 dex at both 
wavelengths. All elements show increasingly more negative 3D-1D corrections with decreasing 
metallicity, larger in magnitude at X = 8500 A. They are the largest in the case of magnesium and 
iron, and are approximately equal to -0.80 dex at [M/H] = -3.0 for X = 8500 A and the excitation 
potential of eV (Fig. 2). 

3 http://www.aip.de/~mst/linf or3D_main.html 
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Figure 2: 3D-1D abundance corrections for neutral Figure 3: Same as in Fig.2 but at A = 8500 A. 
atoms plotted as function of metallicity and excitation 
potential at A = 4000 A. 



The smallest 3D-1D abundance corrections in the absolute sense were observed for the highest 
excitation potentials (Fig. 2 and 3). This is because spectral lines with higher excitation potential 
form in deeper layers of the stellar atmosphere, where differences between the 3D and ID tempera- 
ture stratifications and the amplitude of the horizontal temperature fluctuations are smaller (Fig. 1). 

At solar metallicity the slightly hotter mean temperature of the 3D atmosphere (Fig. 1) reduces 
the line strength of the neutral atoms and horizontal temperature fluctuations increase the line 
strength [8]. These effects make opposite contributions of approximately equal magnitude to the 
line strength in the 3D model, thus producing small 3D-1D abundance corrections. However, 
at low metallicity ([M/H] = -3) the relative importance of horizontal temperature fluctuations is 
significantly larger which produces stronger lines in 3D and thus results in substantial negative 
3D-1D abundance corrections. For neutral atoms, the trends of 3D-1D abundance corrections 
with metallicity are in general agreement with those obtained by [4], despite slight differences in 
the atmospheric parameters of the model atmospheres used in both studies. 

4.2 Abundance corrections for ions 

The dependence of 3D- ID abundance corrections on metallicity for singly ionized atoms is 
shown in Fig. 4 and 5, at A = 4000 A and 8500 A, respectively. The 3D-1D abundance corrections 
are significantly lower than in case of neutral atoms and typically do not exceed ± 0. 1 dex in the 
metallicity range [M/H] = 0.0 to -3.0 (Fig. 4 and Fig. 5). 

Temperature inhomogeneities in the deeper atmosphere layers (log Tr oss w 0) are the main 
cause of negative 3D-1D abundance corrections for the highest excitation potential lines. On the 
other hand, low excitation potential lines are mostly influenced by different temperature stratifica- 
tions predicted by the 3D and ID models in the outer atmospheric layers, although the differences 
between the two are small here (Fig. 1). 

It is important to note that 3D-1D abundance corrections for ionized atoms in our work be- 
come increasingly more negative with increasing excitation potential at fixed metallicity. This is in 
contrast to the findings obtained by [4]. 
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Figure 4: 3D-1D abundance corrections for singly Figure 5: Same as in Fig.4 but at A = 8500 A. 
ionized atoms plotted as function of metallicity and 
excitation potential at X = 4000 A. 



5. Conclusions 

We have investigated the influence of convection on the line formation in the photosphere of 
red giant stars, focusing on the comparison of abundances derived using the 3D hydrodynamical 
(C0 5 BOLD) and ID hydrostatic (LHD) stellar atmosphere models. We have found that the mag- 
nitude and sign of the 3D- ID abundance corrections strongly depends on stellar metallicity and 
atomic parameters of a given spectral line, such as wavelength and excitation potential. At the 
lowest metallicities ([M/H] < -2.0) the differences between the predictions of the 3D and ID mod- 
els may become as large as -0.80 dex for neutral atoms and -0.15 dex for singly ionized atoms. 
The large difference in abundances predicted by the 3D and ID models at low metallicities are 
caused by temporal and horizontal temperature inhomogeneities in the 3D models and different 
temperature stratification predicted by the 3D and ID stellar atmosphere models. Impact of the 
temperature inhomogeneities and lower temperatures in the outer layers of 3D model atmospheres 
of metal-poor giants in the context of spectral line formation was first noted by [1], later discussed 
by [4] and confirmed here using a different 3D code. We therefore conclude that if the predictions 
of currently available 3D stellar atmosphere models are indeed correct these results may signal a 
warning regarding the usage of ID stationary stellar atmosphere models in stellar abundance work 
at low metallicities ([M/H] < -2.0). 

The 3D-1D abundance corrections obtained in this work both for neutral and ionized atoms 
are in general agreement with those derived by [4]. However, the abundance corrections obtained 
by us are smaller, especially at the lowest metallicity. We argue that significantly larger temperature 
difference in the outer atmospheric layers between the 3D and ID model atmospheres used by [4] 
compared with the models used in this work may be partly responsible for this discrepancy. 

In view of these somewhat discrepant findings, we would like to point out that the overall 
methodology of the two studies is somewhat different since our 3D and ID models are constructed 
using the same micro-physics (opacities, equation of state) which is not the case in [4]. 
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